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Abstract In this paper, we first discuss tlie detection of supernova neutrino on Eartli. Tfien we propose a 
possible metfiod to acquire information about S13 smaller than 1.5° by detecting the ratio of the event numbers 
of different flavor supernova neutrinos. Such an sensitivity cannot yet be achieved by the Daya Bay reactor 
neutrino experiment. 
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1 Introduction 

The supernova (SN) explosion is one of the most 
spectacular cosmic events and a source of new physics 
ideas'^'-'— •— Observable efTects of SN neutrinos 
in underground detectors have been a subject of in- 
tense investigation in astroparticle physics, both on 
general grounds and in relation to the SN event like 
1987A. In particular, flavor oscillation in the SN may 
shed light on the problem of neutrino masses and mix- 
ing by means of the associated matter effects. Several 
neutrino laboratories, including the Daya Bay reac- 
tor neutrino underground laboratory's^ which is un- 
der construction, can be used to detect possible neu- 
trino events from an SN explosion and serve as the 
SN Earth Warning System. Hence theoretical pre- 
diction for the detection of SN neutrinos in the Daya 
Bay and other neutrino experiments is very desirable. 

In the realistic case, when we detect neutrinos 
from a type II SN explosion at Daya Bay, there 
are three effects (the collective effects arising from 
neutrino-neutrino interactions'^' -' -' -' -' iS- Ji- i2J ^ 
the well-known Mikheyev-Smirnov-Wolfenstein 
(MSW) effectsli^'— •— and Earth matter 
effects'^'— need to be considered. Using the 
Landau-Zener formula'^'^, the expression of the 
crossing probability, Ph, which is the neutrino 
jump probability from mass eigenstate Vi to 1^3 at 
the high resonance region inside the SN, can be 



calculated'Sa.^'^. With the relation between Ph 
and the mixing angle ^13, we can predict the SN neu- 
trino event numbers as a function of ^13. Therefore, 
we can propose a possible method to acquire infor- 
mation about small 6*13 through the detection of SN 
neutrinos'iSJ. 

2 Detection of SN neutrinos on Earth 

SN are extremely powerful explosions in the uni- 
verse which terminate the life of some stars'^'— '—'^ . 
They make the catastrophic end of stars more mas- 
sive than 8 solar masses leaving behind compact rem- 
nants such as neutron stars or black holes which may 
be observed. For historical reasons, SN are divided 
into two wide categories (type I and type II) charac- 
terized by the absence or presence of hydrogen lines. 
However, the most important physical characteristics 
is the mechanism that produces the SN, which distin- 
guishes SN of type la from SN of type lb, Ic, and II. 
This difference becomes noticeable in the light spec- 
trum some months after maximum luminosity, when 
the ejecta become optically thin and the innermost 
regions become visible: the spectrum of SN la is dom- 
inated by Fe emission lines, while SN lb, Ic, and II 
show O and C emission lines. From the point of view 
of neutrino physics, type lb, Ic, and II SN arc much 
more important than type la SN, simply because they 
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produce a huge flux of neutrinos and antineutrinos of 
all flavors. 

The type II SN is thought to be generated by 
the core collapse of red (or blue as SN1987A) gi- 
ant star with a mass between about 8-9 and 40- 
60 solar masses. The total energy release (about 
3 X 10^^) is approximately the gravitational binding 
energy of the core. It generates intensive neutri- 
nos which take away about 99% of this total en- 
ergy. The explosion itself consumes about 1% of 
this total energy. The vast amount of neutrinos are 
produced in two bursts. In the first burst which 
lasts for only a few milliseconds, electron neutri- 
nos are generated via the electron capture by nuclei 
e~+N{Z,A)^ N{Z-l,A) + iy^ and the inverse beta- 
decay e~-fp— > n+v^. In the second burst which lasts 
longer, neutrinos of all flavors and with a being 
e, r) are produced through the electron-positron 
pair annihilation e~ + e+ ^ i^^t + i^a, electron- nucleon 
bremsstrahlung + N ^ + N + I'a + i^a, nucleon- 
nucleon bremsstrahlung N + N ^ N + N + + i^a, 
plasmon decay 7 — > i^a + i^a and photoannihilation 
7-1-e* e^+Vc+vJ^^^. When SN neutrinos of a 
definite flavor are produced they are approximately in 
an effective mass eigenstate due to the extremely high 
matter density environment. While they propagate 
outward to the surface of the SN they could experi- 
ence collective effects'^'-'-'-'-'-'— and MSW 
effects^i^'— '— '-i^. After travelling the cosmic dis- 
tance to reach Earth, the arriving neutrinos are mass 
eigenstates, which then oscillate in flavors while go- 
ing through Earth matter. Therefore, we also have to 
consider Earth matter effects'^'— when we com- 
pute the event numbers of the various flavors of neu- 
trions. 

Let P^^ represent the collective effects of neutrino- 
neutrino interactions which is a stepwise flavor con- 
version probability of neutrinos at a critical energy 
Ec ■ In order to obtain a simple expression for P^^ for 
the neutrino and P^^ for the antineutrino, we take a 
constant matter density and box-spectra for both the 
neutrino and antineutrino'^'—'^. An analysis of the 
collective effects in the case of three flavors has been 
made in Ref. 11| , and it allows us to characterize the 



collective oscillation effects and to write down the fla- 
vor spectra of the neutrino and antineutrino arriving 
at Earth. Following Ref. [llj, we have P^^ — P^„ — 1 
for the normal hierarchy; and P^^ = 1, while 



P. 



(1) 



Jl {E<Ec), 
(0 {E>Ec), 

for the inverted hierarchy, where Ec — 7MeV^'^ . 



There are two MSW resonance regions: the high 
resonance region and the low resonance region. Let us 
denote the probability that the neutrinos jump from 
one mass eigenstate to another at the high (low) res- 
onance layer by Ph (Pl)- Using the Landau-Zener 
formula'^'^, one can obtain that 



exp(-^7P)-exphf7(-%r)] 



l-exp[-f7(-%r)] 



(2) 



7 = 



I Am? J sin^26li 



2E cos2e'i3 \dliiNjdr\ 



where is the electron density and F can be calcu- 
lated by Landau's method'^. Using the SN matter 
density proflle p f^i C ■ {10^ cm / r)^ ■ 10^" g / cm^ where C 
is a structure constant between 1 and 15 12^' and 
Psa 1 (small ^13), we can obtain an simple expression 
oiP„ that'3^'^'^'^ 



Ph = exp < 




(3) 



where jArngjI = 2.6 x 10~^eV^. Similarity, we can 
calculate the expression of the crossing probability at 
the low resonance region inside the SN, P^. However, 
due to the large angle solution of the neutrino mixing, 
Pl~0. 




Fig. 1. Illustration of the path of the SN neu- 
trino reaching the detector underground in 
Earth. D is the location of the detector, 6 
is the incident angle of the neutrino, O is the 
center of Earth, L is the distance the neutrino 
travels through Earth, and x is the distance of 
the neutrino to the center of Earth. 

Suppose a neutrino reaches the detector at the in- 
cident angle 9 (see Fig. 1). Then the distance that 
neutrino travels through Earth is 



L — {—RE + h)cos( 



{Re -hy sin'' 



(4) 
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where h (« OAkm for the Daya Bay experiment) is 
the depth of the detector in and Re — 6400 km is 
the radius of Earth. Let x be the distance that the 
neutrino travels into Earth, then the distance of the 
neutrino to the center of Earth, i, is given by 

X = ^/{-Re + hf + (L - xY + 2{Re- h){L - x) cos9. 

Let Pie be the probabihty that a neutrino mass eigen- 
state Vi enters the surface of Earth and arrives at the 
detector as an electron neutrino i^., one obtains'^' ^ 



da;F(x)sinC-..,, (5) 



P2e = sin" 6*12 + i sin' 26*12 / ^ 

where 612 = 32.5°, the potential V{x) that 
the electron neutrino experiences in Earth is 
V^G E pix) / (mp + m„) and is defined as 



dxArr,{x) 



A„(a;) = 



Ar 



2E 



(cos 26'i2-£(2:))2+sin2 26*12, 



where e{x) = 2i?F(a;)/Am2i and p{x) is the matter 
mass density inside Earthly . 

In the following, we calculate the event numbers 
N{i) of SN neutrinos that can be observed through 
various reaction channels " i" . This is done by inte- 
grating over the neutrino energy E of the product 
of the target number Nt, the cross section of each 
channel cr(i), and the neutrino flux function at the 
detector F^{E)/AnD^, 



N{i) = NT / dE-a{i) 



1 



■f! 



(6) 



where a stands for the neutrino or antineutrino of a 
given flavor, and D (10 kpc in the present discussion) 
is the distance between the SN and Earth. After a 
straightforward calculation, the fluxes at the detector 
can be obtained'^' ^ (a; = /i,r): 

pDiM) ^ P2ePHFW + (l-P2eP^)F(f), 

F,T' = (l-P2e)Fi:"+P.eFi!', 

= (l-P2ePH)Pi!> + (l+P2ePH)Pi!\ 

2FmN^ ^ P2ePi!' + (2-P2e)Pi!', (7) 

for the normal hierarchy (Arngj > 0), and 

pn,n = fP2.Pil^' + (l-P2.)Pi°', iE<Ec) 

iE>Ea) 

F,y^ - P^(l-P2e)Pi!^ + (l + P2.P.i-P.^)Pil", 
' (1 - P2.)P(«) + (1 + P2.)P(1') , iE<Ec) 

2F,y^ = (l + P2eP«-P^)P<f' 

+ (1 + P^-P2,P«)P(°\ (8) 



for the inverted hierarchy (Amgi < 0). In Eqs. ([7]) 
and ((HI), Pj^' is the time- integrated neutrino energy 
spectrum of flavor a in vacuum which can be de- 
scribed by the Fermi-Dirac distribution 



T{0) 

a \ J TP rj 



P„3Tiexp(£;/T,-77j + l' 
where is the temperature of the neutrino'^'— 

T,^ = 3 - 4MeV, P^.. = 5 - 6MeV, 



(9) 



T,^ = 7~9MeV, 



(10) 



is the pinching parameter of the spectra to repre- 
sent the deviation from being exactly thermal!^' ^ 

r/..«3-5, 77^^ « 2.0 -2.5, 77,^ « 0-2, (11) 

P^"' is the the luminosity, and F^j is defined by 



P„ 



exp(a;-r/„) + l 



dx. 



In the next section, using the relation between the 
event number of SN neutrinos detected at Daya Bay, 

and the mixing angle ^13, we will propose a possi- 
ble method to acquire information about ^13 smaller 
than 1.5°. 

3 Acquire information about ^13 
smaller than 1.5° at Daya Bay 

It can be seen from the above section that, using 
Eqs. (d]), (P, dSD, ®, O, and we can obtain the 
relation between the event number of different flavor 
SN neutrinos, N, and the mixing angle 6*13. Let R 
be the ratio of the event number of 1/^ over that of 

one can obtain P as a function of ^13. Therefore, 
we can propose a possible method to acquire some 
information about the mixing angle ^13 by detecting 
SN neutrinos. In this section, we will just consider 
the process of the neutrino-carbon reactions, since 
our method is not suitable for the inverse beta-decay 
and the neutrino-electron reactions in the Daya Bay 
experiment. The details are discussed at length in 
Ref. 0. 

The Daya Bay Collaboration uses LAB as the 
main part of the liquid scintillator and the total de- 
tector mass is about 300 tons. LAB has a chemical 
composition including C and H . In our calculation, 
the ratio of the numbers of C and H, Nc/Njj, is 
about 0.6. Therefore, the total numbers of target 
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protons A^;f^ —1.32x10^^. For the neutrino-neutrino 
reactions, the effective cross sections are as follows'—': 

{a{^^C{y,,e-f^N)) = 1.85 x lO'^^ cm^, 
(a(^'C(^>,,e+)''B)) = 1.87xlO-^'cm^ (12) 

for the charged-current capture, and 

= 1.33x10-^^ cm^ 
((T(j>f C)) = 6.88 X 10-^^cm^ (13) 
(cr(j^,(/>,)'^C)) =3.73x lO-^'cm^, x^^i.t, 

for the neutral-current capture. Using Eqs. ([T]), ([3]), 
®, dl]), O, dHl), (HH) and (US]), we plot the ratio R 
of the event number of u^, over that of which can 
be detected at Daya Bay as a function of the mixing 
angle ^13. The result is shown in Fig. 2. 



e = 30° e - 90° 




01 2 3 1 2 3 



( in units of degrees ) 9^^ 



e-93° 9-150" 




1 2 3 1 2 3 



8,3 "13 

Fig. 2. The ratio of the event number of to 
that of S'e, R, as a function of the mixing angle 
013 in the channel of neutrino-carbon reactions 
at the Daya Bay experiment. The incident an- 
gle is (a) e = 30°; (b) 6 = 90°; (c) 9 = 93°; 
(d) 9= 150°. The solid curves correspond to 
the normal hierarchy (max) , the dashed curves 
correspond to the inverted hierarchy (max), 
the dot-dashed curves correspond to the nor- 
mal hierarchy (min), the dotted curves corre- 
spond to the inverted hierarchy (min), where 
"max" ("min") corresponds to the maximum 
(minimum) values of Ta and rja . 



It can be seen from Fig. 2 that the uncertainties 
of R due to Ta and rja are not large. For ^13 < 1.5°, 
R is very sensitive to On- However, while 6*13 > 1.5°, 
R is nearly independent of ^13. Therefore, when ^13 
is smaller than 1.5°, we may restrict the mixing an- 
gle 9-13 in a- small range and get information about 
mass hierarchy by detecting the ratio of event num- 
bers of SN neutrinos even though there are still some 
uncertainties due to the incident angle 9, the mass hi- 
erarchy ArTigj^ , and the structure coefficient C of the 
SN density function. 

At the Daya Bay experiment, the sensitivity of 
sin^ 26*13 will reach 0.01, i.e., to determine 6*13 down 
to about 3°. Therefore, if the actual value of 6*13 is 
smaller than 3°, the Daya Bay experiment can only 
provide an upper limit for ^13. However, if an SN 
explosion takes place during the operation of Daya 
Bay, roughly within the cosmic distance considered 
here, it is possible to reach a much smaller value of 
6*13 through the ratio of the event numbers of differ- 
ent flavor SN neutrinos in the channel of neutrino- 
carbon reactions as discussed above. It is interesting 
to note that because of the multi detectors set up, 
experiments such as the Daya Bay have an internal 
coincidence check for SN neutrino events. 

4 Summary and discussion 

In this paper, we first discuss the detection of SN 
neutrinos on Earth. Since neutrino flavor conversions 
inside the SN depend on the neutrino mixing angle 
^13, we give a possible method to acquire information 
about 013 smaller than 1.5° by detecting SN neutrinos 
at Daya Bay. 

We let the parameters in the neutrino energy spec- 
tra (the temperatures and the pinching parameters) 
vary in some reasonable ranges. In fact, the simu- 
lations from the two leading groups, the Livermore 
groupl^^' and the Garching group'^'^, led to pa- 
rameters which agree within about 20-30%. However, 
their central values of SN parameters are different. 



References 

1 Kotake K, Sato K, Takahashi K, Rcpt. Prog. Phys., 2006, 
69: 971—1143 

2 Dighe A S, Physics potential of future supernova neutrino 
observation, arXiv:hep-ph/0809.2977, Neutrinos from a 
core collapse supernova, arXiv:hep-ph/0712.4386 

3 [Daya Bay] GUO X H et al., yl precision measurement of 



the neutrino mixing angle 6\z using reactor antineutrinos 
at Daya Bay, arXiv:hep-ex/0701029 

4 DUAN H Y, Fuller G M, QIAN Y Z, Phys. Rev. D, 2006, 
74: 123004 

5 DUAN H Y, Fuller G M, Carlson J et al. Phys. Rev. D, 

2006, 74: 105014 

6 DUAN H Y, Fuller G M, Carlson J et al. Phys. Rev. D, 

2007, 75: 125005 



No. X 



GUO Xin-Heng et al: Supernova Neutrinos Detection On Earth 



5 



7 DUAN H Y, Fuller G M, Carlson J, Comp. Scic. & Disc, 
2008, 1; 015007 

8 Hannestad S, Raffelt G G, Sigl G et al. Phys. Rev. D, 2006, 
74: 105010 

9 Raffelt G G, Smirnov A Y, Phys. Rev. D, 2007, 76: 081301 

10 Raffelt G G, Smirnov A Y, Phys. Rev. D, 2007, 76: 125008 

11 Dusgupta B, Dighe A, Phys. Rev. D, 2008, 77: 113002 

12 Chakraborty S, Choubey S, Dasgupta B et al. J. Cosmol. 
Astropart. Phys., 2008, 09: 013 

13 Wolfenstein L, Phys. Rev. D, 1978, 17: 2369—2374 

14 Wolfenstein L, Phys. Rev. D, 1979, 20: 2634—2375 

15 Mikheyev S P, Smirnov A Y, Sov. J. Nucl. Phys., 1985, 42: 
913—917 

16 Mikheyev S P, Smirnov A Y, Nuovo Cimcnto C, 1986, 9: 
17—26 

17 Dighe A S, Smirnov A Y, Phys. Rev. D, 2000, 62: 033007 

18 GUO X H, YOUNG B L, Phys. Rev. D, 2006, 73: 093003 

19 GUO X H, HUANG M Y, YOUNG B L, Realistic Earth 
matter effects and a method to measure small S13 in the 
detection of supernova neutrinos, arXiv:0806.2720 [hep-ph] 

20 Zener C, Proc. R. Soc. London A, 1932, 137: 696—702 

21 Landau L D, Lifshitz E M. Quantum Mechanics: Non- 
relativistic Theory. New York: Pergamon, 1977. 199 — 237 

22 Kuo T K, Pantaleone J, Rev. Mod. Phys., 1989, 61: 937— 
979 

23 Lunardini C, Smirnov A Y, J. Cosmol. Astropart. Phys., 
2003, 06: 009 



24 Kachelrie/3 M, Strumia A, Tomas R et al. Phys. Rev. D, 
2002, 65: 073016 

25 Giunti C, Kim C W. Fundamentals of Neutrino Physics and 
Astrophysics. New York: Oxford, 2007. 511—539 

26 Shapiro S L, Teukolsky S A. Black Holes, White Dwarfs, 
and Neutron Stars. New York: Wiley, 1983. 512—542 

27 Mohapatra R N, Pal P B. Massive Neutrino in Physics and 
Astrophysics, 3nd ed. Singapore: World Scientific, 2004. 
340—358 

28 Brown G E, Bethe H A, Baym G, Nucl. Phys. A, 1982, 
375: 481—532 

29 Bethe H A, Rev. Mod. Phys., 1990, 62: 801—866 

30 Janka H T, Astron. Astrophys., 2001, 368: 527—560 

31 Dziewonski A M, Anderson D L, Phys. Earth Planetary 
Interiors, 1981, 25: 297—356 

32 Janka H T, Hillobrandt W, Astron. Astrophys., 1989, 224: 
49—56 

33 Janka H T, Hillebrandt W, Astron. Astrophys. Suppl. Ser., 
1989, 78: 375—397 

34 Cadonati L, Calaprice F P, Chen M C, Astropart. Phys., 
2002 16, 361-372 

35 Totani T, Sato K, Dalhed H E et al. Astrophys. J., 1998, 
496: 216—225 

36 Keil M T, Raffelt G G, Janka H T, Astrophys. J., 2003, 
590: 971—991 

37 Keil M T, Supernova neutrino spectra and applications to 
flavor oscillations, arXiv: astro-ph/0308228"| 



